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ARTICLE INFO ABSTRACT 
Article history: Under some author's opinion the wind is not a periodical phenomenon and therefore it is more 
Received 22 July 2010 reasonable to invest in renewable periodical energies as tides. In this paper we have developed a 


Accepted 9 September 2010 computer application based in MatLab®, that through the FFT (Fast Fourier Transform) analyzes the 


variations of wind speed amplitude in the time and frequency domain. The data were sampled every 
Keywords: 10 min in the period 2000-2008. The data come from 31 Automatic Meteorological Stations (EMAs), the 


a 4 d country of Mexico and correspond one per state. The survey shows the representation of spectral- 
als a temporal surfaces to long time intervals, as one year or more and denotes seasonal envelopes that 


alter the pattern at certain times everyday. As a conclusion, the wind has an important periodical 
component for the country of Mexico, since the fundamental component of the wind speed 
represents a frequency of 1/24 h~! in a very accurate form throughout the time studied. To harness 
the wind potential of the country of Mexico it should be kept in mind that there is a minimum wind 
speed between 8 and 16h and a maximum close to 24h. 
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1. Introduction quality of life, have developed significant growth in renewable 
power generation across the world [1]. The importance of 
The necessity of mitigating adverse environmental effects implementing tools to support policies to promote the dissemi- 
from the use of fossil fuels, volatile fuel prices and improve the nation of these technologies is now a universally accepted 
standard [2]. The wind energy is one of the most powerful and 
profitable way to satisfy the demands of sustainable development 
* Corresponding author. [3], while the production of wind energy has grown 30% per year 

E-mail address: fmanzano@ual.es (F. Manzano-Agugliaro). over the last 10 years [4]. 
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The wind is considered as one of intermittent renewable energy 
[5], and for this reason energy storages systems are needed [6], in 
terms of energy storage, being the production of hydrogen one of 
the more efficient [7]. 

Denny [8] explains that the wind has a variable and 
unpredictable nature which results in increased challenges for 
electricity system operators. 

Ma et al. [9], say that the random character of wind is 
significant since the wind speed depends on the terrain and the 
height, and the unpredictability and variability of the wind power 
generation is one of the fundamental difficulties faced by power 
system operators. 

Wind studies to evaluate the capacity of the wind potential of a 
geographic area are of two types: physical models and statistical 
models. Physical methods have advantages in long-term predic- 
tion whereas statistical methods do it well in short-term 
prediction [9]. 

Examples of statistic models are applied in Argelia [10], in 
Turkey [11] or in Minnesota, USA [12]. In Mexico there has been 
a recent work to determine the wind potential of the state of 
Veracruz [13] and to predict the wind in the state of Oaxaca [14]. 

The technique of the spectral analysis through the Fast Fourier 
Transform (FFT) allows to determine whether a phenomenon has 
periodic components, based on calculating its harmonic [15]. The 
use of the FFT representing spectral-temporal surfaces has been 
used in oceanography, in particularly in the study of hurricane 
Wilma in Mexico, to determine the temporal evolution of waves, in 
particular variations associated to the period and trajectory of the 
waves [16]. 

The purpose of this paper is to analyze, by applying the FFT 
representation of spectral-temporal surface, the amplitude varia- 
tions of wind speed in the time and frequency domain, in order to 
determine whether there exists a periodic behavior. 


2. Materials and methods 
2.1. Data 


In this work we have used the data from the automatics 
meteorological stations (EMAs) in the country of Mexico [17]. We 
have analyzed 31 stations, one per state, the geographic distribu- 
tion is shown in Fig. 1 and names in Table 1. The data used are a 
series of records from January 1st 2000 to December 31st 2008, 
with an interval of 10 min between them. 

In Fig. 2 it is shown an example of the data used, corresponding 
to Presa Emilio Lopez Zamora meteorological station (Baja 
California Norte state) from 1st-01-2007 to 31st-12-2007. 


2.2. Methods 


2.2.1. Spectral analysis of the wind 

An effective technique for studying the periodicity of natural 
processes is the analysis in the frequency domain, which 
determines the distribution of energy depending on the frequency 
or speed of repeated procedures. When the information is discrete, 
as the one obtained through sampling techniques that is the most 
common method, treatment is usually based on the use of discrete 
Fourier transform or FFT in its optimized version for computing 
speed. Thus, the Fourier spectrum H(f) of the register denoted by 
h(t) was obtained by the formula [18] 


H(f) = I j h(t) eet 


where fis frequency, t is time andj is the imaginary unit (j = V—1). 
In general, H(f) is a complex quantity called Fourier transform of 
h(t). H(f) can also be expressed as 


Kilometer 
3600 


Fig. 1. Geographic distribution of EMA. 
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Table 1 
FT Analysis of EMAs in Mexico. 
ID Name station 
Cal Calvillo 
PELZ Presa Emilio Lopez Zamora 
CSL Cabo San Lucas 
CDC Ciudad del Carmen 
Pal Palenque 
Bas Basaseachi 
NRo Nueva Rosita 
ENCB Escuela Nacional de Ciencias Biológicas 
Ame Agustín Melgar 
PAII Presa Allende 
Cal Ciudad Altamirano 
Pac Pachuca 
Tiz Tizapan 
PMa Presa Madín 
Ang Angamacutiro 
IMTA Instituto Mexicano de Tecnología del Agua 
Aca Acaponeta 
PEC Presa El Cuchillo 
PAn Puerto Ángel 
UTT Universidad Tecnológica de Tecamachalco 
Hui Huimilpan 
Can Cancún 
CVa Ciudad Valles 
Obi Obispo 
Her Hermosillo 
Par Paraíso 
Alt Altamira 
Hua Huamantla 
Jal Jalapa 
Cel Celestún 
Zac Zacatecas 


State Period obtained (Hz) 
Aguascalientes 1.157e—005 
Baja California Norte 1.157e—005 
Baja California Sur 1.157e—005 
Campeche 1.157e—005 
Chiapas 1.157e—005 
Chihuahua 1.157e—005 
Coahuila 1.157e—005 
Distrito Federal 1.157e—005 
Durango 1.157e—005 
Guanajuato 1.157e—005 
Guerrero 1.157e—005 
Hidalgo 1.157e—005 
Jalisco 1.157e—005 
México 1.157e—005 
Michoacan 1.157e—005 
Morelos 1.157e—005 
Nayarit 1.157e—005 
Nuevo León 1.157e—005 
Oaxaca 1.157e—005 
Puebla 1.157e—005 
Querétaro 1.157e—005 
Quintana Roo 1.157e—005 
San Luis Potosi 1.157e—005 
Sinaloa 1.157e—005 
Sonora 1.157e—005 
Tabasco 1.157e—005 
Tamaulipas 1.157e—005 
Tlaxcala 1.157e—005 
Veracruz 1.157e—005 
Yucatan 1.157e—005 
Zacatecas 1.157e—005 


H(f) = R(f) + jF) = |H( fy? 


where R(f) is the real part of H(f) and I(f) is the imaginary part. The 
amplitude spectrum of h(t) is denoted by |H(f)| and the phase 
spectrum is denoted by ®(f). These are defined as 


IH(f)|l = y R( f) + P(A) D(f) = tan" GP) 


R( f) 


When the amplitude spectrum of h(t) is studied, it is 
possible to identify the principal harmonic components of the 


Wind Speed (m:s-!) 


| Mij | 


digital record, because they always appear with major ampli- 
tude [16]. Frequencies corresponding to principal harmonic 
components are the frequencies of the tide waves present in the 
register. These waves are identified visually and their period T 
(T=1/f) can be read directly from the amplitude spectrum. 
Thus, these harmonic components are easily identified in 
the amplitude spectrum. Obviously, spectral analysis can be 
used to identify and quantify any harmonic component in a 
digital register. For this reason, the spectral analysis is 
considered today a standard tool of analysis in many scientific 
fields. 


Wii iy 
ý i aN 1 


Year 2007 Data 


Fig. 2. Example of the data (Presa Emilio Lopez Zamora in Baja California Norte state from 1st-01-2007 to 31st-12-2007). 
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Tm 


Window width k samples, interval [0,Tm-k] s 


k samples FFT 


espec [t=Tm-K-0] 


FFT 


espec [t=Tm-k-1] 


Fig. 3. Outline of sw developed to calculate the FFT. 


2.2.2. Characterization of the average day 

To characterize the average day for each station in the years 
2007 and 2008, we have developed an Excel software application 
that calculates the average of each 10 min interval for the 24 h to 
finally obtain the average day of the year. This information is 
important for the statistical comparison of areas. 


3. Theory/calculation 


The study of the frequency of wind requires the capture and 
management of a large amount of data, a sampling period of a few 
minutes and a period of observation of the order of years to 
appreciate daily and seasonal periodicity. 

The direct application of FFT to the entire data block captured 
from the wind will not provide very valuable information. It will be 
necessary to develop a specific process of treatment of these 
signals that denote the characteristic parameters that can be 
expected, such as if the existence of periodicity and what type of 
variations occur along a period of a year. 


Amplitude 


Sl . 
o 1 Fred 


A sampling period Tm of 10 min in the capture of data allows to 
collect the most important information of the variable, but 
generates sets of the order of 50 KS per year. 

The direct application of the FFT to the set of samples would 
provide spectral information that would not reflect seasonal 
changes, so we have developed the following specific technique of 
analysis. 

The set of samples collected of the variable is fragmented by 
means of windows of time of appropriate type. A rectangular 
window introduces distortions due to convolution with Tm-Sinc(f 
Tm) in the domain of frequency. Distortion increases the 
importance of a less number of samples of the window [19]. 
Window Hamming is the most suitable for this application. The 
temporal width of the window is set to: Tm-k, obtaining a set of 
unidimensional vectors of size k of the variable. Then we apply the 
FFT to each vector, which produces two vectors of size k, 
corresponding to the module phase respectively. By means of a 
process of decimation that eliminates the elements of the FFT 
image we get two vectors of size K/2: espec[t=tm-k-O], 
espec[t = tm-k-1] as shown in Fig. 3. 


Fig. 4. Sample of spectral-temporal surface for Presa Emilio Lopez Zamora. 
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Fig. 5. Sample of FFT analysis for Presa Emilio López Zamora. 


After that we generate a two dimensional array Z through the 
expansion of the vectors by rows spec [t=tm-K-0], spec 
[t = tm-K-1]..., corresponding to the FFT module: 


es pec(t = tm-k-0) 
Z- es pec(t =tm-k-1) 


es pec(t =tm-k- j) 


Each row vector of the array Z corresponds to the time interval 
Tm-k, 2 Tm-k,... and thus, for each element zij, j and i correspond to 
the time and frequency variables, respectively. The array Z allows a 
representation by a surface as shown in Fig. 4, which corresponds 
to an interpolated surface built by the Matlab application, using the 
surfl function. In one of the axis, time is represented, the second 
one corresponds to the frequency and amplitude is represented in 
the axis z. 


This type of treatment allows the study of the spectrum for a 
given interval of time and the visualization of the set for the study 
of the behaviour for a long time or seasonal. The spectral projection 
of a given period reflects the periodicity of the peaks that are 


produced using the same frequency value regardless of the period 
and place of the measure. 


4. Results and discussion 
4.1. Spectral analysis 


Applying the methodology described above we get a spectral- 
temporal surface for each station from the variations of amplitude 
of the wind speed in the time and frequency domains. In Fig. 4 we 
show a typical example of these surfaces for the Emilo Lopez 
Zamora Reservoir meteorological station, Nuevo Leon state, where 
there is a concentration of energy at the lowest frequencies. This 
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Fig. 6. Average day for the North area of Mexico. 
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Fig. 7. Average day for the Pacific area of Mexico. 


indicates that the phenomenon is not random a priori. Then, to 
prove if there exist periodicity, it is studied a spectral-temporal 
surface section in the frequency domain. 

In Fig. 5 it is shown an example at the same meteorological 
station, where the analysis through FFT shows a peak above the 
noise in the frequency of 1.157e—005 Hz, which is found in all the 
EMAs studied in the country of Mexico, see Table 1. This frequency 
corresponds to a period of 24 h. The other peaks are consequence of 
the design of the software made to show that a function is not 
sinusoidal. 

The survey of the wind periodicity has required capture and 
manages of a large amount of data, with a sampling period of 
10 min and an observation time of order of years to be able to 
observe daily periodicities. 

This implies that an envelope seasonal changes the daily 
pattern in certain periods. In these isolated cases we have detected 
temporal intervals of high frequency variability of wind speed that 
are currently under study, although abnormalities were discarded 
because the EMA have been presented in different EMAs. 
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4.2. Average day 


Once obtained an average day in every season, these are 
represented in four graphics, grouped by areas within the country 
of Mexico: North (Fig. 6), Pacific (Fig. 7), center (Fig. 8) and Gulf of 
Mexico (Fig. 9). 

It is observed that in the north zone, Fig. 6, the average day of 
the year 2007 has a behavior similar to that one in 2008, i.e., the 
curve has approximately the same shape. It is also observed that in 
all these stations, the wind speed has a minimum value between 8 
and 16h, and a peak around the 24h, which shows the possible 
phase shifting of the eolic energy generation with the possible 
daily energetic necessities. 

On the other hand it can also observed that the average day of a 
year is higher than the average day of other year depending on the 
season, and there does not exist a year with more speed wind than 
others for every season in this area. We detach the identic 
behaviour of the average day for the station Nueva Rosita during 
the two analyzed years. This may implie that the geographic 
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Fig. 8. Average day for the Center area of Mexico. 
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components of this station located in the Durango desert are not so 
important in the variation of the speed wind year over year than in 
the rest of the stations. 

For the Pacific area, Fig. 7, it is shown a similar behaviour to the 
North area in the minimum values of the speed wind, excepting for 
the station of Cape San Lucas, probably due to its special location in 
the peninsula of Baja California, where it is shown only a minimum 
at 16:00. Here all the stations have higher speed wind than in 2007, 
excepting the mentioned of Cape San Lucas. On the other hand the 
station of Hermosillo, located in the Sonora desert presents an 
identic behaviour in the two analyzed years, such as it was shown 
the station of Nueva Rosita. 

In the Center area, Fig. 8, the speed wind presents the same 
minimum from 8:00 to 16:00 than in the other two mentioned 
areas. Here there does not exist a year with higher speed wind in all 
the stations. However, we detach that there exist less variations in 
the speed wind year over year than in the other two areas studied 
till the present, being the speed wind more stable in this area. Here 
we also have an station that shows an identical behaviour during 
the two years, Huamantla, that it is not in a desert but it has a 
common properties with Hermosillo and Nueva Rosita, and this is 
that it has a flat geography. 

In the Gulf of Mexico area, Fig. 9, there also exist a minimum for 
8:00 to 16:00 in the wind speed, excepting Ciudad del Carmen 
station, whose behaviour is stable during the whole day. Paraiso 
and Cancun stations show a behaviour year over year less similar 
between each other that the rest of stations in this area. We detach 
the difference between the speed wind during the years 2007 and 
2008 for Ciudad del Carmen station, showing both of them the 
same behaviour due probably to the geographic situation, i.e., both 
are located in an island. 


5. Conclusions 


In this survey it has been demonstrated that the methodology 
used through the Fast Fourier Transformed (FFT) in the time and 
frequency domain is useful to analyze the possible periodic 
components of the wind speed. The harmonics component obtained 
through the FFT, besides the fundamental, indicate that the function 
associated to the wind speed respect time is not sinusoidal. 

The survey shows the representation of spectral-temporal 
surfaces to long time intervals, as one year or more, denote 


Hour Hour 


e Gulf of Mexico area. 


seasonal envelope that alter the pattern at certain times 
everyday. 

The wind have important periodical components in the country 
of Mexico, such as the fundamental component of the wind speed 
that presents a frequency of 1/24h7!, very accurate along the 
period of time studied. 

To harness the wind potential of the country of Mexico it should 
be kept in mind that there is a minimum wind speed between 8 
and 16h and a maximum close to 24h. 


Acknowledgements 


The authors gratefully acknowledge the contributions of 
Emmanuel Alvarez Ramirez to provide the data of Automatic 
Meteorological Station in Mexico. This work has been funded by 
Project TEC2009-13763-C02-02 (Ministry of Science and Innova- 
tion of Spain). 


References 


1] Cook NJ. Towards better estimation of extreme winds. J Wind Eng Ind 

Aerodynamics 1982;9:295-323. 

2] Sebitosi AB, Pillay P. Renewable energy and the environment in South Africa: a 

way forward. Energy Policy 2008;36:3312-6. 

3] Sebitosi AB. Energy efficiency, security of supply and the environment in South 

Africa: moving beyond the strategy documents. Energy Policy 2008;33: 

1591-6. 

4] Pacala S, Socolow R. Stabilization wedges: solving the climate problem for the 

next 50 years with current technologies. Science 2004;305:968-72. 

5] Moriarty P, Honnery D. Intermittent renewable energy: the only source of 

hydrogen? Int J Hydrogen Energy 2007;32:1616-24. 

6] Barton JP, Infield DG. Energy storage and its use with intermittent renewable 

energy. IEEE Trans Energy Convers 2004;19:441-8. 

7] Hall PJ, Bain EJ. Energy-storage technologies and electricity generation. Energy 

Policy 2008;36:4352-5. 

Denny E. The economics of tidal energy. Energy Policy 2009;37:1914-24. 

Ma L, Luan SY, Jiang CW, Liu H, Zhang Y. A review on the forecasting of 

wind speed and generated power. Renew Sustain Energy Rev 2009;13: 

915-20. 

[10] Himri Y, Rehman S, Draoui B, Himri S. Wind power potential assessment for 
three locations in Algeria. Renew Sustain Energy Rev 2008;12:2495-504. 

[11] Ucar A, Balo F. A seasonal analysis of wind turbine characteristics and 
wind power potential in Manisa, Turkey. Int J Green Energy 2008;5: 
466-79. 

[12] Wichser C, Klink K. Low wind speed turbines and wind power potential in 
Minnesota, USA. Renew Energy 2008;33:1749-58. 

[13] Cancino-Solorzano Y, Xiberta-Bernat J. Statistical analysis of wind power in the 

region of Veracruz (Mexico). Renew Energy 2009;34:1628-34. 


oO œ 


728 Q. Hernandez-Escobedo et al. / Renewable and Sustainable Energy Reviews 15 (2011) 721-728 


[14] Cadenas E, Rivera W. Short term wind speed forecasting in La Venta, 
Oaxaca, México, using artificial neural networks. Renew Energy 2009;34: 
274-278. 

[15] Brigham EO. The fast Fourier transform and its applications. New Jersey: 
Prentice HallNew Jersey; 1988. 

[16] Silva-Casarin R, Mendoza-Baldwin E, Escalante-Mancera E, Marino-Tapia I, 
Ruiz-Renteria F. Wind waves induced by hurricane Wilma in Puerto 


Morelos, Quintana Roo, Mexico. Ingenieria Hidraulica en Mexico 2009;24: 
93-109. 

[17] SMN, Servicio Meteorológico Nacional, México. Visited 14th July 2009, http:// 
smn.cna.gob.mx/productos/emas/. 

[18] Bath M. Spectral analysis in geophysics. Amsterdam: Elsevier; 1982. 

[19] Capuano S, Harris FJ. An introduction to FFT and time domain windows. 
Instrum Measur Mag IEEE 2007;10:32-44. 


